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Abstract We propose the use of quantum polarization spectroscopy for detecting 
multi-particle entanglement of ultracold atoms in optical lattices. This method, 
based on a light-matter interface employing the quantum Farady effect, allows 
for the non destructive measurement of spin-spin correlations. We apply it to the 
specific example of a one dimensional spin chain and reconstruct its phase diagram 
using the light signal, readily measurable in experiments with ultracold atoms. 
Interestingly, the same technique can be extended to detect quantum many-body 
entanglement in such systems. 



1 Introduction 

Accurate measurements of quantum correlations in the next generation of exper- 
iments with ultracold atoms in optical lattices are one of the most challenging 
obstacles in the quest for the realization of quantum simulators of magnetic sys- 
tems '' 3 . Not only local order parameters, but also long range correlations are nec- 
essary for the faithful discrimination of magnetic phases in the strongly correlated 
regime. One of the advantages of using optical lattices for simulating solid state 
systems, is that atoms are easily coupled to light with extremely accurate control. 
The atomic matter properties are then inferred from the measurement of the light 
scattered off of an atomic sample. 

Here, we review a proposal for inferring spin-spin correlations using a quan- 
tum polarization spectroscopy scheme based on a light matter interface 4 . The idea, 
put forward initially in Ref. 5 and further developed in Refs. 6,7 , consists in cou- 
pling the polarization of a beam illuminating the optical lattice with the magnetic 
momenta of the trapped atoms. As we explain in this work, using a light beam in a 



1 :Ffsica Teorica: Informacio i Processos Qualities, Universitat Autdnoma de Barcelona, E-08 193 
Bellaterra, Spain 

2:ICREA, Institucio Catalana de Recerca i Estudis Avancats, E08011 Barcelona 



2 




Fig. 1 (Color online) Schematic detection setup, atoms placed in an optical lattice of periodicity 
d/2 (thin line; red lattice) are illuminated by a laser beam in a standing wave configuration (dark 
line; blue lattice) shifted by a from the optical lattice configuration. The output light is redirected 
through a polarimeter which measures its polarization through a homodyne detection (HD). 



standing wave configuration, as sketched in Fig. Q] one can achieve a modulation 
of the light-atoms coupling which allows one to reconstruct the atomic spin-spin 
correlations. One of the most important advantages of this method is that it is non 
destructive, i.e. the atomic sample is kept in the trap and can be reused for further 
measurements. We apply this method to a specific one-dimensional spin chain. We 
show how to measure the model phase diagram by accessing the order parameters 
of the different phases. 

Furthermore, we will show that apart from measuring spin-spin correlations, 
polarization spectroscopy allows to discriminate whether a magnetic system, in 
our case a spin chain, is entangled or not. This is a long standing problem in 
the context of quantum information theory and many-body systems (see for ex- 
ample 8 ), since it is extremely difficult to characterize quantum entanglement for 
general many particle states. In some specific cases, such as spin systems, one can 
derive spin squeezing inequalities involving the system total angular momentum 
which reveals whether the many-body state is entangled or not. In this context, 
several proposals have been put forward based on particle scattering, e.g. neu- 
tron scattering, to probe entanglement in magnetic systems 9 ' 10 ' 11 . Here we show 
that entanglement witnesses based on spin-squeezing inequalities are straightfor- 
wardly measured with our proposed scheme. Moreover, the flexibility of opti- 
cal setups in modulating the periodicity of the probe wave give a lot of freedom 
and accuracy in the measurements compared to neutron scattering. We expect this 
scheme to open a new route for the, so far elusive, detection of many-body entan- 
glement. 

The paper is organized as follows: in Sec.|2]we review the quantum polariza- 
tion spectroscopy technique aimed at measuring spin-spin correlations in optical 
lattices while in Sec. [3] we discuss how to simulate spin chains in optical lattices; 
in Sec.|4]we show the numerical results for the output signal for the reconstruction 
of the phase diagram of the spin model; in Sec.[5]the entanglement detection using 
quantum polarization spectroscopy is described, and finally in Sec.|6]we conclude. 
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2 The detection scheme 

The detection scheme based on polarization spectroscopy that we use in this work 
has been described in^&i. The scheme consists in shining the atoms with a non 
resonant probe beam in a standing wave configuration as shown in Fig. [T| Due 
to the Faraday effect, the polarization of the incoming light is rotated as a conse- 
quence of an effective magnetic field generated by the atomic magnetic moments. 
By measuring the change in the polarization of the output light we acquire infor- 
mation on the total angular momentum of the atomic sample. 

For light propagating along the z-axis, parallel to the atomic array, light-atom 
interaction is best expressed using the Stokes parameters defined as: 

s i = ^(a'la x - a'l,a y ) , (1) 

*2 = \^a\,a x + a\a y ) , (2) 

S3 = ^-X a l a x-ala y ), (3) 



where a x and a y are the photon annihilation operators with polarization along x 
and y such that n x = a\a x and n y = ala y are the number of photons per unit of 
time with polarization x and y respectively. Using this definition, the atom-light 
interaction is described by the Hamiltonian: 

H AL = ~Ks 3 J e / f , (4) 

where the coupling constant jc depends on the optical depth of the atomic sample 
and on the probability of exciting an atom due to the probe. The effective angular 
momentum depends on the intensity profile of the probe beam. In the case of 
a simple standing wave it is given by 

/j = /= y, c nS Z n , (5) 
V-L n 

and the coefficients are defined as 



c n = lj dzcos 2 [kp(z-a)]\w(z-nd)\ 2 , 



(6) 



where kp is the wavevector of the probe light, a is a shift and w(z — nd) is the 
first band Wannier function of the atom centered at lattice position z = nd. In 
the calculations, for simplicity, we approximate the Wannier functions with delta 
functions centered at the lattice positions so that the coefficients are now given by 
c n = 2cos 2 [kpd(n — a)} where we defined the dimensionless shift a = a/d. 

We assume the incoming light to be strongly polarized along the x direction, 
i.e. (Si) = N p h S> 1 where Si = J dtSj and N p h is the beam total number of photons. 
Therefore we can approximate the other two Stokes operators as two effective 
conjugated variables: X = S2/ y/Nph ar, d P = S3 / \/Nph sucn mat 

[X,P) = § L ~i- (7) 
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Integrating out the Heisenberg equations of motion for these light quadratures we 
obtain 

X out =X in -Kjff, (8) 

where Xj n is the quadrature of the incoming light, and X out is the output light 
emerging from the sample that can be measured using homodyne detection as 
shown in Fig. [TJ Since we assumed the initial beam to be strongly polarized along 
the x direction, (Xi„) = 0, we obtain: 

(X 0Ut ) = -K{J e /f), (9) 

thus, the mean of the effective angular momentum is mapped into the mean output 
light quadrature. Similarly, higher moments (fluctuations) are also mapped. In this 
way all the moments of can be extracted by the noise distribution of the output 
light and in particular the variance: 

(AX ollt ) 2 = (AX in f + K 2 {Ajff)\ (10) 

where (4X„,) 2 is the input noise (for a coherent state (AXi n ) 2 = 1/2). Differ- 
ent magnetic phases can be distinguished by studying the mean effective angu- 
lar momentum {Jt ) an d the variance (AJ Z ^^) 2 . The former immediately tells 
us whether the spin chain is ferromagnetic or not. The second gives us access to 
magnetic correlations: 

e(k P ,a) = [AJ e / f f = 

= — ^cos 2 [k P d(m - a)]cos 2 [k P d(n - a)} 

nm 

x %{m,n). (11) 

where & z (m,n) = (S zm S zn ) — (S zm )(S zn ) is the two-point correlation function. As 
noticed in 6 , in the case of a sample for which the net magnetization is zero, the 
output signal can be connected to the magnetic structure factor. Indeed by averag- 
ing the signal over the phase shift a one gets: 

e{k P ) = J dae{k P ,a) = -S(2k P ) (12) 

In principle, averaging over the phase shift a we are losing some information 
on the correlations. In this work we assume an accurate control on the shift a and 
we define the quantity: 

Ae(kp,a h a 2 ) = e(k P ,ai) -e(k P ,a 2 ) (13) 

which is the difference of the signal with fixed wavevector kp at two different 
phase shifts. We will show that the quantity Ae(kp, (X\ , Ctq), by appropriately choos- 
ing the parameters kp,0Ci,CC2, can be linked to the local order parameters necessary 
for the identification of the different phases of a given model. 

If we assume no net magnetization, then the expression for Ae simplifies to: 

Ae(k P ,al,a2) = -!-^{cos[2fcpt/(m + n-2ai)] 

- cos[2kpd(m + n-2a 2 )}}& z {m,n) (14) 



5 



3 Realization of spin-1 Hamiltonians in optical lattices 

Here we briefly review how to simulate spin chains with ultracold atoms in opti- 
cal lattices and discuss its phase diagram. Spin-1 atoms confined in a deep opti- 
cal lattice generated by two counter-propagating lasers of wavelength X are well 
described, within the tight binding approximation, by the Bose-Hubbard Hamil- 
tonian 13 . Defining the creation and annihilation operators aj a and a, jCT in site i 
of an atom with spin components a = 1,0,-1 along the quantization axis, the 
Hamiltonian takes the form: 

Hbh = ^ £«/(«/ + 2m) - ii £«; + (15) 



~ f L( a lCT a <+'. (7+/j - c -) 



The operator n, = Y.a a ]a a >,o ls the total number operator of site i while S; — 
Y, a a'^in^a a' a i& i s the spin operator (matrices T are the usual spin-1 angular 
momentum operators and we use h = 1). The parameters appearing in the Hamil- 
tonian Hbh are: the usual Hubbard repulsion Uq, a spin dependent interaction U2, 
the chemical potential /J and the tunneling rate t. While the chemical potential 
fixes the total number of atoms, the remaining parameters can be evaluated from 
the depth of the optical lattice and from the scattering lengths associated with 
different scattering channels—. 

The phase diagram of this model in the fi~t plane consists of insulating lobes 
as in the spinless Bose-Hubbard model where the lobes size depends on the ratio 
£/2/£/o— • For ur, it filling and for sufficiently small tunneling t the system is in a 
Mott insulator state with one atom per site. Virtual tunneling of the atoms between 
neighboring sites gives rise to an effective magnetic interaction described by the 
bilinear-biquadratic Hamiltonian—: 

H BB = /£cos(e)S; • S i+1 + sin(0)(5 / • S i+1 ) 2 (16) 
i 

The Hamiltonian ( 116b is derived within second order perturbation theory in the 
ratio t/U a , Of = 0, 2 and the relevant parameters read: 

^ A- (17) 

2 



where the angle 6 varies in the interval [— k; n}. 

Hamiltonian dl6l) is characterized by a rich phase diagram, sketched in Fig. [2] 
depending on the angle 9 and which has been extensively studied in the liter- 
ature, se c t 6 ! 17 ! 24 ! 18 ! 19 - 20 ! 15 - 21 ! 22 ] 23 and references therein. Here we briefly discuss 
the model phase diagram and the corresponding order parameters. 

The ferromagnetic phase.- For n/2 < 6 < 5k/4 the ground state is ferromag- 
netic: all the spins, breaking the rotational symmetry of Hbb, align along some 
direction with a net spontaneous magnetization, which serves as a local order 



6 




Fig. 2 Phase diagram of the bilinear-biquadratic Hamiltonian i l!6t in the interval 8 6 [— lt;n\. 
The four phases are: the ferromagnetic phase (F), the critical phase (C), the Haldane phase (H) 
and the dimer phase (D). 



parameter. For the remaining values of 9 the ground state lacks of spontaneous 
magnetization. However, within this interval we can distinguish different phases. 

The critical phase.- In the interval k/4 < 9 < k/2 the system is in a critical 
phase in which the model is gapless due to soft collective modes at momenta 
q = 0, ±2n/(3d) where d = A/2 is the distance between two adjacent sites. The 
ground state organizes in slightly correlated clusters of three neighboring spins 
(trimers). This fact is reflected in the spin-spin correlation functions (S Z jS z ( i+r \) 

which show aperiod-3 oscillations 24 . In momentum space this feature emerges as 
a peak at q = 2k /(3d) in the magnetic structure factor defined as: 

S(q) = We i ^ m ~"\S zm S m ). (19) 

mn 

Recently Lauchli et al. 21 have shown that nematic (i.e. quadrupolar) correlations at 
momentum q = 2k /(3d) are enhanced in the critical phase while spin correlations 
become smaller when increasing 9 from 0.2k to 0.57T. Together with the absence 
of the gap, the enhanced nematic correlations are a distinctive feature of the critical 
phase. 

The Haldane phase. - For — k/4 < 9 < k/4 the system is in the Haldane phase 
which is gapped and contains for 0=0 the spin- 1 isotropic Heisenberg chain and 
for tan(0) = 1/3 the Affleck-Kennedy-Lieb-Tasaki (AKLT) point for which the 
ground state is exactly known 16 . Numerical results in this region based on den- 
sity matrix renormalization group (DMRG) simulations show that decreasing 9 
from k/4 to 9l ~ 0.13147T, the so called Lifshitz point, the peak at momentum 
q — 2k / (3d) in the magnetic structure factor moves continuously to q = n/d (see 
Ref. 19 ). Notice that, although these peaks signal some correlations, the presence 
of a gap excludes local long range magnetic order and spin correlations decay ex- 
ponentially. The Haldane phase can be instead characterized in terms of a hidden 
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topological order parameter, called the string order parameter' 

% (m,n) = ( S zm exp in V S zJ S zn ) (20) 




This order, being topological, cannot be revealed with local measurements. 

The dimer phase.- The interval — 37T/4 < < — n/4 is still debated. At 6 — 
—n/4 the gap closes and for smaller values of 6 it reopens again. In this region the 
ground state breaks translational invariance and organizes in slightly correlated 
dimers. For — 3n/4 < 8 < — n/2 it is still under debate whether the system is 
always dimerized or it becomes nematic as proposed by Chubukov 17 . Numerical 
result s 18 i 20 i 15 i 21 show that the dimer order parameter: 

D=\(Hi-H i+1 }\ (21) 

where //, = cos(0)S, • 5,+ i +sin(0)(S, • S,+ i) 2 , is different from zero up to values 
very close to 6 — —3n/4 giving strong evidence for the absence of the nematic 
phase except only in an infinitesimal region close to 9 = —3n/4. 



4 Phase diagram reconstruction 

In this section we discuss the results of the detection scheme applied to the bilinear- 
biquadratic Hamiltonian. The quantities e(kp,a), which depend on all possible 
correlations between two spins, are computed numerically by means of the DMRG 
algorithm 12 . We simulate spin chains with open boundary conditions and lengths 
which are multiple of 2 and 3 reducing known finite size effects due to incom- 
mensurability—. In the DMRG simulations we choose the number of block states 
sufficiently large to ensure that the truncation error is less than 10 . 

The ferromagnetic phases is easily detected by looking at the average value of 
the effective angular momentum: 

(Jff(k = 0)) = ^='£s zn (22) 

which is proportional to the total magnetization along the z direction. 

Since (Jf ) is zero in the other three phases, we need the second moment of 
J e z ^ in order to characterize this phase. In Fig. [3] we show e(kp, a) in the critical, 
Haldane and dimer phases. A common feature of the three phases is the presence 
of a high peak at kpd = n/2 due to antiferromagnetic correlations. Apart from this, 
the three plots are qualitative different. In fact for 9 > 8l, the Lifshitz point, the 
signal is characterized by peaks at kpd ~ tt/3 and kpd ~ 2n/3. These resemble 
the peaks of the magnetic structure factoiQ and are due to the period-3 oscillations 
of the correlation functions. We will study these correlations in Sec. 14. H and show 
that they detect the critical phase. For 9 < —n/4 we find the appearance of other 
small peaks at kpd = n/4 and kpd = 3n/4 signaling a different order with a larger 
period. We will study more carefully these features in Sec. 14.21 

1 Notice that from Eq. (Tf2j, e(kp,a) is related to the structure factor S(2kp) at double the 
value of the momentum. 
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Fig. 3 Left column, the function e[kp, a) for different values of 9 in the three phases for L = 
132: top 6 = — 0.57T (dimer), middle 0=0 (Haldane), bottom 9 = 03k (critical). Right column, 
the same plots but restricted to a = 0. 



Since the presence of these distinctive peaks is relevant for the determination 
of the phase of the spin chain, we find it convenient to subtract the background 
generated by all possible correlations in definition (lilt by instead using the quan- 
tity Ae(k Pl d\ , (X2) defined in Eq. (TT3l . 

To see how to choose the parameters kp, 0C\ , OL2, let us consider the dimer phase. 
In this case we find it convenient to choose kp = n/Ad which is the periodicity of 
the dimers. Then we study the behavior of e(jc/Ad, a) in one point of the dimer 
phase as a function of a as shown in Fig. [4] The quantity e(n/Ad, a) is an oscil- 
lating function of a. In order to optimize the information on the correlations at 
kp = n/Ad we choose the difference between the maximum at a.\ =3/2 and the 
minimum at Oq. = 1/2. Thus as an indicator of the critical phase, we will study 
the quantity Ae(n/Ad, 3/2, 1 /2). In the critical phase, a similar analysis leads to 
k P = n/Ad, ai =5 /A, a 2 = l/2 (see alsol). 



4. 1 Detecting the critical phase 

For the critical phase, we have seen that the distinctive peaks are at kpd = n/3, 
and as shown in the previous section we choose (X\ =5/4 and Ct2 = 1 /2. Thus we 
define the quantity 



C E =4e(w/3d,5/4,l/2) = 
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Fig. 4 The quantity s(n/4d, a) for 6 = —0.5n (dimer phase) for L = 132 as a function of a. 



where we used the fact that the ground state is an eigenstate of the total angular 
momentum with zero eigenvalue: 



The quantity C £ is sensitive to correlations which oscillate with a period 3 and 
represents a footprint of the critical phase. In fact, in Fig. [5] we show the signal 
C £ for different values of 9 in the antiferromagnetic phase between —O.ln and 
0.57T. The results clearly show that the critical phase is very well detected by a 
positive value of C £ . For 9 = 0.2n, in the Haldane phase and close to the phase 
transition, we still observe a large positive value, probably due to residual period 3 
correlations persisting in the Haldane phase for 9 > 9t. However for 9 = 0.2n we 
find a non negligible dependence with the size of the sample. A finite size scaling 
analysis suggests that in the thermodynamical limit for L — > °° the quantity C £ goes 
to zero as 1 /L for 9 = 0.27T, while for the other values of 9 > 0.24 n it converges 
to a finite value (see Ref. 7 ). Our findings indicate that by measuring C £ which 
depends only on spin-spin interactions we are able to infer the occurrence of the 
phase transition and thus the quantity C £ behaves as an order parameter for the 
critical phase. 



4.2 Detecting the dimerized phase 

Let us now consider the dimerized phase. As discussed before the presence of 
peaks at kp = n/4d signals pairing of neighboring spins. Notice that if we average 
the signal e(kp, a) over a these peaks disappear. Therefore these features are not 
visible in the magnetic structure factor. 




(23) 



^j(S zm S zn ) — {S zm y.'SV. 



(24) 



n n 
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Fig. 5 (Color online) The quantities C e = Ae(7l/3d,5/4,l/2) (squares) and D £ = 
Ae{n/Ad, 1/2,3/2) (circles) as a function of G for L = 132. We distinguish the model phases 
with different shading: horizontal lines (dimer), no shading (Haldane), oblique lines (critical). 
The solid and dashed lines are only guides to the eye. 



We find that the quantity 



D £ = Ae{n/4d, 1/2,3/2) 



41 



sin 



71 ' I 

12 {M 



& z (m,n) 



(25) 



is suitable for the detection of the dimer phase. The factor sin [n/2(m + «)] ensures 
that only the pairs of spins with positions m and n of opposite parity contribute to 
D £ . Moreover the sin function gives an alternating sign depending on whether the 
distance between the sites is even or odd. Therefore the quantity D E is an extension 
to long range correlations of the dimer order parameter D defined in Eq. (121b . 

In Fig.[5]we show the results for the signal D £ for different values of 0. Similar 
to the dimer order parameter D, the quantity D E is significantly different from 
zero only in the dimerized phase, therefore acting as an alternative dimer order 
parameter. 



5 Entanglement detection 

Detecting entanglement in many-body systems is not an easy task. In magnetic 
systems, such as the spin chain considered in this work, one can employ spin 
squeezing inequalities based on collective angular momentum operators (see>& for 
a review). An entanglement witness is an operator which is positive valued for all 
separable (non entangled) states, while there exists at least one entangled state for 
which the expectation value of the witness is negative. The witness we propose is 
based on the effective angular momentum defined in ©. The construction follows 
Refs.M 
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As before we define an effective angular momentum which we assume can be 
written on the light fluctuations: 

J a = ^c m S am a=x,y,z (26) 
in 

where now we consider the angular momentum fluctuations in the two other di- 
rections. Let us consider the quantity: 

V= £ AJ 2 a = £ ZciCjdSaiSa^-iSadiSaj)) (27) 

a=x,y,z a=x,y,z ij 

Now if the many-body system is in a product state: 

Pprod=P\®Pl®---®pN (28) 

we have: 

(SaiSaj) - (S ai )(S aj ) = 8ij ((S 2 ai ) - (S ai ) 2 ) (29) 
Using the relation for spin 5 particles: 

(4> + (4> + (4)=^ + l) (30) 

and the inequality: 

(S xi ) 2 +(S yi ) 2 + (S zi ) 2 <s 2 (31) 
we see that for product states: 

Vprod^S^C 2 (32) 

i 

If we consider separable states: 

Psep = Y J PnPn,sep, < p n < 1, £/>h = 1 (33) 

n n 

where each state p n>sep in the mixture is separable, we have 

V sep = £44 > £>„£(44)„ > £p„,£ c 2 = s^cf (34) 

a n a n i i 

where the first inequality comes for a mixture p = Y^nPnPn- AX 2 > p n (AX 2 ) n 
and (AX 2 ) n is the variance evaluated in the nth ensemble element; the second 
inequality comes from Eq. (132b . 

Therefore, a possible entanglement witness is given by the quantity: 

W=y-j£c? (35) 

i 

Notice that the coefficients c, and consequently the quantity V depend on the probe 
light momentum k and on the shift, a, between the optical lattice and the probe 
light on the standing wave configuration. Both parameters can be changed giving, 
therefore, an important and necessary flexibility for detection of different entan- 
gled states. In Fig. [6] we show W for a = as a function of k for states in the 
critical, Haldane and dimer, phases. It is evident that for certain values of k, W is 
negative thus detecting entanglement. 

This method provides an operational entanglement detection scheme which is 
scalable, robust and that can be realized in present-day experiments with ultracold 
atoms in optical lattices. We stress that the quantity W is very general and can be 
used even if the sample is subject to thermal fluctuations or disorder. 
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Fig. 6 (Color Online) Expectation value of the entanglement witness W from Eq. 05t for three 
values of 6 in the three different phases: G = —0.5k (solid (red), dimer phase) G = 0.102^: 
(dashed (green), AKLT point in the Haldane phase), G = 0.3k (dotted(blue), critical phase). In 
the numerical simulations we take L = 96. All these states are clearly detected for small enough 
values of k. 



6 Conclusions 



We have presented a probing technique based on matter-light interface for the in- 
vestigation of quantum correlations in magnetic systems simulated by ultracold 
atoms in optical lattices. We have shown that this scheme permits to obtain exper- 
imentally the order parameters of non trivial magnetic phases by homodyne mea- 
suring the fluctuations of the probing light quadratures after crossing the atomic 
sample. Moreover, we have shown that this technique, which is not destructive and 
realizable with present technology, allows also to detect experimentally the entan- 
glement in magnetic non trivial many body systems without actually carrying out 
an unnecessary state tomography. 
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